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Abstract: We perform a detailed analysis of baryon asymmetry generation in the z/MSM 
(an extension of the Standard Model by three singlet Majorana fermions with masses below 
the Fermi scale) . Fixing a number of parameters of the i/MSM by the neutrino oscillation 
data, we determine the remaining domain of the parameter space from the requirement of 
successful baryogenesis. We derive, in particular, the constraints on the mass splitting of a 
pair of singlet fermions, and on the strength of their coupling to ordinary leptons, essential 
for searches of these particles in rare decays of mesons and in beam-dump experiments 
with intensive proton beams. 
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1. Introduction 

There are many ways the baryon asymmetry of the Universe (BAU) could have arisen. They 
differ from each other by the nature of baryon and lepton number non-conservation, CP- 
violation and by the manner a departure from thermal equilibrium is realised. For practical 
purposes it is important to distinguish between the scenarios which can or cannot be 
experimentally tested. By experimentally testable we mean those in which particle physics 
experiments (limiting the time-scale by the existing accelerators including the LHC, and by 
the ILC, CLIC or an 67 collider for the future) could potentially determine all parameters 
of the theory, necessary for theoretical computation of the amplitude and of the sign of the 
baryon asymmetry of the Universe^. 

The number of experimentally testable scenarios is quite limited. Basically, all the 
mechanisms for baryogenesis, which use new physics at energies much higher than the 
Fermi scale do not fall in this category. For example, the GUT baryogenesis (for a review 
see Q) is not testable, as it includes superheavy particles with masses of the order of 
0(10^^) GeV, the direct search of which is impossible. The thermal leptogenesis (for a 
review see 0) is only partially testable. It did pass the important tests coming from the 
observed neutrino masses and mixings, but the prediction of the baryon asymmetry in it 

^Of course, the requirement of testability is not a "physical" one in a sense that the Nature is not obliged 
to be nice to us. 
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depends on CP-breaking phases which cannot be measured in low energy experiments, as 
the typical mass of the superheavy Majorana leptons is 0{W^^) GeV - too high to allow 
their creation in laboratory. 

If new particles, responsible for BAU, are in the reach of existing or future colliders, 
the baryogenesis mechanism has more chances to be testable. One of the examples is the 
resonant leptogenesis of [^, ^. It uses essentially the same physics as thermal leptogenesis, 
but relies on a possible degeneracy of heavy Majorana leptons. This allows to shift their 
masses to the region accessible by the LHC or linear colliders, but still larger than the 
electroweak scale. Though the interactions of one of the Majorana leptons, relevant for 
BAU, are too weak to lead to its discovery at colliders two other heavy leptons can be 
found and a peculiar flavour structure of their interactions can provide the "smoking gun" 
signature of the resonant leptogenesis § ^ . 

Yet another example is the electroweak baryogenesis (for reviews see ||7|, ^, |9|) where 
the source of baryon number non-conservation is the rapid high-temperature anomalous 



processes with fermion number violation |10|, and departures from thermal equilibrium 



are due to strongly first order electroweak phase transition [T^, The electroweak 

baryogenesis requires an extension of the Standard Model right above the electroweak scale, 
leading to specific predictions for LHC and ILC (see [^, |l^ for recent discussions). 

The baryogenesis can also occur due to new particles with the masses considerably 
smaller than the Fermi scale. An example (actually, the the only one, known to the 
authors) is associated with the minimal extension of the Standard Model (SM) in neutrino 



sector - the z^MSM |15, 16|. In this model one adds to the SM three singlet Majorana leptons 
with masses smaller than the electroweak scale and uses the standard sea-saw Lagrangian, 
though with the completely different choice parameters (for a review see p!7[| ) . The lightest 
neutral singlet lepton with the mass in 0(10) keV region can be the dark matter (DM) 
candidate; the other two with masses M ~ 0{1) GeV would then generate the baryon 
asymmetry of the Universe. The Higgs boson with non-minimal coupling to gravity in the 
same model can play the role of the inflaton and make the Universe flat, homogeneous 
and isotropic, producing cosmological perturbations leading to structure formation |18|. 
The fact that so modest and natural modification of the SM can address at once many 
problems which the SM cannot solve (neutrino masses and oscillations, dark matter, baryon 
asymmetry of the Universe and inflation) forces us to take it seriously and study in detail 
the predictions of z^MSM, which can be used as a guideline for experimental searches for 
new particles it contains^. 

The astrophysical observations coming from X-ray satellites tell that the coupling of 
the lightest (DM) singlet fermion A'^i to active neutrinos is so small that it practically does 



not contribute to the see-saw formula [24|. This means that two other singlet fermions 



of the i/MSM A''2,3 must be able to explain simultaneously the pattern of neutrino masses 



^We thank A. Pilaftsis for discussion of this point. 

^Further extensions of the uMSM are surely possible. See |20| |2^] for discussion of the model 
with extra scalar singlet and for analysis of the model with higher-dimensional operators added. Of 
course, the predictions of these models are less certain than those of the j^MSM. 
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and oscillations and lead to baryogenesis. The parameter counting goes as follows^. Out 
of 11 new parameters of the model 7 can be (potentially) measured in experiments with 
active neutrino flavours only (these are 3 mixing angles, one Dirac CP-violating phase, 
one Majorana phase, and 2 active neutrino masses, the lightest one is massless in this 
approximation). The other 4 can be conveniently chosen as M = ^(Mi + M2) - the 
average Majorana mass, AMm = ^(Mi - M2) - the mass splitting, e < 1 - the ratio of 
the strengths of the couplings of N2 versus N3 to leptons, and finally r/ ~ an extra CP- 
violating Majorana phase. Whether N2^3 can be found experimentally depends crucially 



on 2 parameters - their mass and e |25]. The main aim of the present paper is to find the 
constraints on these two parameters from the requirement of successful baryogenesis. 

The idea that a pair of light (with masses in GeV region) and almost degenerate 
Majorana leptons can result in baryogenesis through their oscillations goes back to ref. p6|| . 
The equations governing baryogenesis in the z^MSM, incorporating the relevant physical 
effects (singlet fermion oscillations, transfer of fermion number from singlet species to active 
leptons and back, etc.) were formulated in |l^ , where it was shown that the i^MSM can lead 
to baryon asymmetry of the Universe and give a dark matter candidate, being perfectly 
consistent with the data on neutrino oscillations. In that paper a perturbative solution 
to the kinetic equations was found, valid in a certain part of the parameter space^. A 
semi-quantitative analytic consideration of baryogenesis in the i/MSM beyond perturbation 
theory has been carried out in |2^, where different essential time scales were analysed, and 
the generic dependence of the baryon asymmetry on some of the z^MSM parameters was 
elucidated. 

An accurate prediction for properties of singlet fermions call, however, for a complete 
exploration of the parameter space and thus for a numerical solution of the kinetic equa- 
tions, without any assumptions about hierarchies of different time scales, made in analytical 
computations. To identify the parameter-space allowing for baryogenesis, we calculate the 
baryon asymmetry as a function of the constants of the z/MSM as follows. We fix the 
parameters of the neutrino mass matrix (^12, O23, ^13, Am^^^, and Am^^^) to their ex- 
perimental values. For any given mass M, mass difference AMm and e we determine the 
maximum of BAU as a function of unknown CP-violating phases and require that this value 
is larger than the observed one. This allows to fix the 3-dimensional domain of parameters 
(M, AMm, e)) most relevant for experimental searches of A^2,3- 

This paper is organised as follows. In Section § we summarise the basic properties of 
the model and formulate the kinetic equations for baryogenesis. In Section ^ we describe 
the results of numerical analysis. The Section ^ is conclusions. A number of appendixes 
contain details of the numerical procedure. In fact, a direct numerical solution of the kinetic 
equations poses a technical problem, as the corresponding differential equations belongs 
to the so-called "stiff" type. The reason is the existence of the many vastly different time 
scales, associated with oscillations and with relaxations of distinct types of deviations from 



*From now on we omit the DM singlet fermion from consideration, as irrelevant for active neutrino 
masses and for baryogenesis. 

^It is required that both A^2,3 are out of equilibrium at the sphaleron freezing temperature Tsph and that 
the rate of A^2,3 oscillations exceeds the rate of the Universe expansion at Tsph- 
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thermal equilibrium. In appendix |A| we transform the kinetic equations, governing the 
evolution of the system into a convenient form, which allows to elucidate the hierarchy 
of the relevant time-scales (the rate of the singlet fermion oscillations, coherence loss, 
transfer of leptonic number from sterile to active species, back reaction, leptonic flavour 
non-conservation, etc) and describe an effective way to solve them numerically. In appendix 
1^ we discuss the dependence of baryon asymmetry on the parameters M, AMm and e, 
and its time evolution. In ^ we compute the maximal baryon asymmetry which can be 
created in z^MSM. 

2. Basic equations 

This chapter is a brief review of the baryogenesis in the i/MSM. First, we will describe the 
model, then we will formulate the kinetic equations, and finally we present the results of 
perturbative computation of the asymmetry, to be compared with numerics later. We use 
the same notations as in [p7|]. 

2.1 Lagrangian and parameters 

The part of the Lagrangian of the z^MSM we are interested at can be written as: 

/ -V2 _ . \ ^ 

C = iNid^-i^'Ni + t\—j^L2N2 + v^e^^'/^^sA^s j ^ 

AMm - 

-MiVliVg ^{NIN2 + NIN3) + h.c, (2.1) 

where and N3 are the singlet leptons and Mib AMj\,f their Majorana masses, = eij^j 
is the Higgs doublet and v = 174 GeV its vacuum expectation value, e is the ratio between 
the strengths of the couplings of and N2 (e < 1 by convention) and r/ is a CP-violating 
phase as explained in the Introduction. The constant J-"^ is expressed through active 
neutrino masses rrii as 

^=^E^^- (2-2) 

i 

The 6 out of the 11 new parameters of the model are explicitly written in the Lagrangian 
( |2.1| ) (M, AMm, r] and rrii). The other 5, related to the active neutrino mixing matrix, 
are hidden is L2 and L3 which are combination of Lg, and Lr, the lepton doublets. 
These parameters are 3 mixing angles 612, ^13, 623, one Dirac phase (p and one Majorana 
phase. It will be denoted by a for the normal active neutrinos mass hierarchy and ^ 
for the inverted one, see | |28[ | for convention. The relations between L2,3 and the flavour 
eigenstates were found in [^], we present them here for ^13 = and 6*23 = f , which is 
within the experimental error bars (for a review see [p8| ). 
For the normal hierarchy we have: 

L^-Lt L^ + Lt 

L2 = +ai — -j= h a2Le + 03 — -j= — , 

Lii — Lt Ln + Lt 

L3 = -ai-i^-^ a2Le + as — , (2.3) 
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where the coefficients are 



ai = ie-^(°+<^)sin^cos0i2 , 



02 = le smf3smfc'i2 

03 = cos Q , 



tan Q 




0.4 . 



For the inverted hierarchy the expressions are similar : 

L 



with the coefficients : 



V2 



^'-7=2 



V2 



cos 6*126 ""^ + i sin 6*126+"^ 



cos ^126"^*^ + i sin ^126 



(2.4) 



(2.5) 



(2.6) 



For the analysis of kinetic equations, it is convenient to define the Yukawa coupling 
matrix Fq,/ between the singlet fermion and the lepton doublet in the following way: 



(2.7) 



which corresponds to a basis where the singlet fermion Majorana mass matrix Ai is 
diagonal and equals to: 



M 



M - AM, 



M 







M + AM, 



M 



(2.8) 



2.2 The framework 

A convenient way to treat a system of interacting Majorana fermions and neutrinos is 



to use a density matrix p |£9|, 31 1. Since we consider 2 singlet fermions and 3 active 
neutrinos and their anti-particles, /9 is a 10 x 10 matrix. For baryogenesis the temperatures 
are above the electroweak scale Tw ~ 100 GeV, we also take M <^ Tyy. In this case 
the matrix p can be simplified [16|: it can be described by 3 chemical potentials pa for 
active neutrino species and by two 2x2 density matrices for the sterile neutrinos, one for 
positive helicity states, pN and one for negative helicity states pN- The evolution of the 
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system satisfies the kinetic equations which take into account creation of singlet fermions, 
their oscihations, and generation of lepton asymmetries |jl^. Introducing the CP-even and 
CP-odd deviations from thermal equilibrium: 



5p^ 
5p- 



PN + PN 



PN — PN, 



P 



eg 



(2.9) 



where p^'^ = exp(^ ~ t) '^i^h p the momentum of the sterile neutrino |1£], these kinetic 
equations are : 



.d6p^ 



d5p. 
\ — — 
dt 



[Re/?; 



-{ReTN,5p+} + i-smip(F'^ p.F - F'^ p.F 



+ '-[lmWj;}\5p.] + -^{lrnTN,6p.], 



[ReH'^\5p. 



^{RePTv, 5p-} + i- sin ^(^F^fiF + F^ fiF' 



+ 2i[lmH'j:^\6p+] + {lmTN,Sp+}, 



.dp 

I— = diag 
dt 

T 
16 



^{rL,p} 



+ i^sinip(^F{6p- + 25p+)F^ - F*{5p^ - 25p+)F'^^ 



Here F^j is Yukawa coupling matrix introduced in the previous chapter, 

T T 
Ftv = — sin 99 F^F, TL = —sinip 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



are the rates of creation of singlet and active fermions correspondingly, and the oscillation 
Hamiltonian is given by 



Tjint 



Lf^F + ^ diag [ - MAMm, MAMm] 

o 1 



(2.14) 



The expressions (|2.13 , 2.14| ) were found for temperatures above Ty/ in the symmetric phase 
of electroweak theory, simp ~ 0.02 is related to the ratio of the absorptive to the real part 
of neutrino propagator in the medium | 



The kinetic equations are written in the so-called monochromatic approximation, which 
assumes that the different energy modes of singlet fermions are effectively decouple. This 
is true in our case, as the collisions between A^2,3 can be surely neglected, and the problem 
is linear in terms of the density matrix for the singlet fermions^. 

The equations (2.10,2.11, 2.12| ), supplemented by initial conditions with no lepton 
asymmetry for the active leptons {p = 0) and zero abundance for the singlet fermions 



-p^'^, Sp- = 0) describe the evolution of active and sterile neutrinos from CP- 



symmetric initial state. The choice of the initial condition for the CP-even part of the 



®In general, this is not the case for active neutrinos in core collapse supernova, where one needs to take 
into account neutrino-neutrino interactions, leading to non-linear effects [ pi] , We thank E. Akhmedov 
and A. Smirnov for discussion of this point. 



-6- 



density matrix 6pj^ is motivated by the following facts. First, the z^MSM interactions of 
singlet fermions are so weak that they are out of thermal equilibrium till temperature T+ 
which is close to the electroweak scale (see the estimates of T+ as a function of the pa- 
rameters of the i/MSM in [^]). Therefore, they are not produced at T > T+. Second, in 
the minimal setup of the i^MSM, where inflation occurs due to the Higgs boson [^], the 
initial conditions for the kinetic equations can be found. They coincide exactly with those 
formulated above |23]. Moreover, the dependence of the produced baryon asymmetry on 
the initial condition for the CP-even part of the density matrix is linear in (5/3+. Even if 
the singlet leptons had relatively strong non-i/MSM interactions at high energy scale, their 
initial (for our kinetic equations) abundance will be in the interval —p'^'^ < 5p+ < 0. In 
other words, the produced baryon asymmetry for generic initial conditions for (5/9+ can only 
be smaller than what we found. So, the prediction of the domain of the parameters, were 
baryogenesis in the j/MSM is not possible, and therefore, identification of the experimental 
goals for their search (Figs. ^ - ^ is robust. At the same time, if |(5/?+| < p^'^, the region 
where baryogenesis in the z^MSM is possible will decrease correspondingly. 

An analytic solution to these equations in a particular limit, discussed in the following 
subsection, was found in The qualitative features of this system of kinetic equations 
have been explored in [^] . We will use them to compute numerically the baryon asymmetry 
of the Universe as a function of the i/MSM parameters in Section |3[ Note that the transfer 
of the lepton asymmetry to the baryon asymmetry because of sphalerons is not included 
into eqns. (2.10|,2.11,p.l2), this can be done analytically in a standard way (see, e.g. [pi]). 



2.3 Analytical expression for baryon asymmetry 



The equations ( p. 101 ), ( |2.11| ) and ( p. 12 ) can be solved perturbatively for small Yukawa 



couplings, provided the following assumptions are satisfied p6[ |: 

(i) Singlet fermions are out of thermal equilibrium for all temperatures above the sphaleron 
freeze-out, Vn ■ tw < 1 with Fat ~ sinc/?^^, tw — Mq/T^^, where Mq = 7 ■ 10^'^ 
GeV. 



(ii) The baryon asymmetry creation temperature Tl — (4MAM^mMo)^/'^ |2^, 16] is much 
higher than the electroweak temperature, Tl > Tyy- 

(iii) The mass difference between two singlet fermions is much larger than the mass dif- 
ference between active neutrinos, AMm ^ 0.04 eV for the normal hierarchy and 
AMm ^ 8 ■ 10~^ eV for the inverted hierarchy. 

The baryon to entropy ratio (the observed one is ^ ~ (8.4 — 8.9) • 10^^^) is given by 

'^c^1.10-^Tt{5p.)\t^ , (2.15) 
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where Tr{6p-.)\T^ can be found in and can be conveniently parametrised as 

^CP = %Y1 Im[F„3[FtF]324j . (2.17) 

Due to the conditions ([|) and (P), the first and second brackets respectively are always 
smaller than 1. This formula was used recently to compute the baryon asymmetry of the 
Universe in the i/MSM, where CP-violation is solely originated from the CP violation in 
the mixing matrix of active neutrinos (i.e. for r/ = 0) in [ p5[ . 

Qualitatively, li T]\f ■ tw S> 1, the singlet fermions equilibrate, and dilute the baryon 
asymmetry. If Tl <C Tw, the oscillations of singlet leptons had no time to develop and the 
creation of baryon asymmetry is strongly suppressed In other words, the maximal 

asymmetry is generated when F^r • tjy ~ 1 and Tl ~ Tw simultaneously. 

In general case, the conditions (|)- are not satisfied, and numerical solution of the 
kinetic equations is necessary. The results of this analysis are presented in the next Section. 

3. Domain of the parameters leading to baryon asymmetry 

In this chapter, we determine the 3-dimensional domain of the parameters e, M and AMm 
that can lead to the observed baryon asymmetry. It is found in the following way. The 
active neutrino mass matrix parameters are fixed to be 

vr 



sin 012 = V0.3, 013 = 0, 023 = ^, 



1712 = \l ^^loi ~ ^ i^^^j iTT-s = y ^f^ttm = ^0 meV, (3.1) 



m2 = = J l\mi^^, (3.2) 



for the normal hierarchy and 



for the inverted hierarchy. The variation of the parameters in ( |3.1| , ^.2| ) within experimental 
error bars (shown, for instance, in |2^) does not lead to any significant changes in the 
admitted parameter range. The unknown CP-violating phases are chosen in such a way 
that baryon asymmetry is maximised (see appendix |C|). They are close to the values given 
below, 

r/ ~ ~ ^, a ~ |, C =^ p (3.3) 
which corresponds to the following Yukawa couplings ratio: 

\Fe2? ■■ \Ff,2\^ : |F,2|' ~ 1 : 11 : 11, (3.4) 
for the normal hierarchy and 

\Fe2\^ ■■ \F^2\^ : |F,2|'=^48: 1 : 1, (3.5) 
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for the inverted hierarchy. The second set corresponds to the model I of ||25[] and the first 
set is somewhat between the models II and III, described in |25|. The sphaleron freezing 
temperature is set at T\y = 140 GeV corresponding to the Higgs mass mn = 125 GeV [Q. 
We have also done the computations for Ty\r = 170 GeV which correspond to mn = 200 
GeV, but the results are very close to the case T\y = 140 GeV. 

With this set of parameters the kinetic equations ( p.l0| ), ( 2.11 ) and ( p. 12 ) are solved 
numerically with the methods described in appendix The admissible domain of the 
parameters is determined from the condition ub / s\theoreticai > ^b/ s\observed- For masses 
above ^ 10 GeV the computation should not be trusted, as the kinetic equations are only 
valid in the limit M <^ Ty/- We expect that the admitted region of the singlet fermion 
masses in the z^MSM closes up at M ~ Mw |27|. A number of intermediate results, such 
as dependence of baryon asymmetry on e, M and AM^ , and its time dependence, can be 
found in appendix 

In the left panel of Fig. || we present the projection of the allowed region to the 
e — M plane for both neutrino mass hierarchies. Here, for every e and M the asymmetry 
is extremised with respect to AMm; in the interior region the asymmetry is greater than 
the observed one, the boundary corresponds to ub/ s\theoreticai = it-b / s\observed- The baryon 
asymmetry generation is possible for very light singlet fermions, M ~ 10 MeV for the 
normal hierarchy and M ~ 1 MeV for the inverted one. They correspond to small e ~ 
few X 10"''. 




Figure 1: Values of e - M (left panel) and AMm - M (right panel) that leads to the observed 
baryon asymmetry for the normal hierarchy (blue - solid line) and for the inverted hierarchy (red - 
dashed line). 

The right panel of Fig. || gives a projection of our domain to the AMm — M plane, 
which is derived in the similar way (asymmetry is now extremised with respect to e at 
fixed AMm and M). The admitted mass difference ranges from a fraction of eV to MeV, 
depending on the singlet fermion mass. 

In Fig. |2| we present yet other projections to the AMm — e plane for several singlet 
fermion masses and two different types of the hierarchies. In general, for larger masses the 
amount of the parameter space where baryogenesis happens is larger. Also, the inverted 
hierarchy has more parameter space than the normal one. 
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Figure 2: Values of AMm and e that leads to the observed baryon asymmetry for different singlet 
fermion masses, M = 10, 100 McV, 1, and 10 GeV. The blue (shortest dashed) line corresponds to 
M = 10 MeV, red (short dashed) - to M = 100 MeV, brown (long dashed) - M = 1 GeV and green 
(longest dashed) to M — 10 GeV. Left panel - normal hierarchy, right panel - inverted hierarchy. 



For most important parameters for experimental searches of the singlet fermions are 
their mass and e, which determines the strength of the coupling of A^2,3 to ordinary leptons. 
Usually it is expressed through the mixing angle U'^ between ordinary and singlet leptons 
(we take a sum over mixing angles with all active neutrino flavours, for an exact definition 



see [25|). For small e. 



= ^ • (3.6) 

The region, where baryogenesis is possible in C/^ — M plane is shown in Fig. ^, which is 
the main result of this paper^. It is derived from Fig. |l] with the use of (|3.6D . We also plot 
there the exclusion regions coming from different experiments such as BEBC [^], CHARM 
|3^, and NuTeV |8| and CERN PS191 experiment ||, g (see also discussion of different 



experiments in [^]). For the case of normal hierarchy, only CERN PS191 have significantly 
entered into cosmologically interesting part of the parameter space of the z/MSM, situated 
below the mass of the kaon. If the hierarchy is inverted, there are some constraints even 
for higher masses. The lower constraint on [/^, coming from baryon asymmetry of the 
Universe, is somewhat stronger than the "see-saw" constraint. In Fig. Q we present the 
expected life-time of the singlet fermions in an experimentally interesting region M < 2 
GeV. 

A detailed discussion of possible experiments and signatures of neutral leptons leading 



to BAU generation in the i/MSM can be found in [25|. For the reader convenience, we 
summarise below the main conclusions of this work. 

Several distinct strategies can be used for the experimental search of these particles. 
The first one is related to production (C/^ effect). The singlet fermions participate in 
all the reactions the ordinary neutrinos do with a probability suppressed roughly by a 
factor U"^. Since they are massive, the kinematics of, say, two body decays — ?■ fJ-^N, 



^These figures supersede other similar plots which appeared previously in a number of works pq , [l7 
and in conference proceeding. 
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Figure 3: Constraints on U'^ coming from the baryon asymmetry of the Universe (sohd hnes), from 
the see-saw formula (dotted hne) and from the big bang nucleosynthesis (dotted line). Experimental 
searched regions are in red - dashed lines. Left panel - normal hierarchy, right panel - inverted 
hierarchy. 




Figure 4: Constraints on the lifetime tn coming from the baryon asymmetry of the Universe 
(solid lines), from the see-saw formula (dotted line) and from the big bang nucleosynthesis (dotted 
line). Experimental constraints from PS 191 are shown in red - dashed lines. Left panel - normal 
hierarchy, right panel - inverted hierarchy. 

— )• e^N or three-body decays Kl,s — )• vr^ -|- -)- A''2,3 changes when is replaced 
by an ordinary neutrino. Therefore, the study of kinematics of rare K, D, and B meson 
decays can constrain the strength of the coupling of heavy leptons. The precise study of 
kinematics of rare meson decays is possible in <I> (like KLOE), charm, and B factories, or in 
experiments with kaons where their initial 4-momentum is well known (like NA48, NA62 
or E787 experiments). 

The second strategy is to use the proton beam dump (C/^ effect). As a first step, 
the proton beam heating the fixed target creates K, D or B mesons, which decay and 
produce The second step is a search for decays of in a near detector, looking for 

the processes "nothing" — )• leptons and hadrons. To this end, quite a number of already 
existing or planned neutrino facilities (related, e.g., to CERN SPS, MiniBooNE, MINOS, 
J-PARC or LBNE at Fermilab), complemented by a near dedicated detector can be utilised. 

For the mass interval M < Mk, both strategies can be used. If mx < M < rriD, the 



- 11 - 



search for the missing energy signal, potentially possible at beauty, charm, and r factories, 
is unlikely to gain the necessary statistics and is impossible at hadronic machines like the 
LHC. Thus, the search for decays of neutral fermions is the most effective opportunity. 
The dedicated experiments on the basis of the proton beam NuMI or NuTeV at FNAL, 
SPS at CERN, or J-PARC can touch a very interesting parameter range for M < 1.8 GeV. 
Experiments like NuSOnG |42] and HiResMz^ [43| should be able to enter in a cosmologically 
interesting region for masses and mixing angles of singlet fermions. 

Going above D-meson but still below S-meson thresholds is very hard if not impossible 
with the present or planned proton machines or B-factories. To enter into a cosmologically 
interesting parameter space would require the increase in the present intensity of, say, 
CERN SPS beam by two orders of magnitude or to produce and study the kinematics of 
more than 10^° B -mesons. For M above the beauty threshold ~ 5 GeV the baryogenesis 
in the i/MSM becomes untestable. 



4. Conclusions 

In this paper we found the 3-dimensional domain of z^MSM parameters (M, AMm, e) that 
may lead to observed baryon asymmetry of the Universe. It ranges from the singlet fermion 
masses in MeV region to tens of GeV and from mass differences from a fraction of eV to 
MeV. Thus, baryogenesis is a generic consequence of the i^MSM with nearly degenerate 
Majorana fermions. Our results can be used to determine an ultimate goal for searches 
for cosmologically interesting heavy neutral leptons, responsible for neutrino masses and 
oscillations, and for baryon asymmetry of the Universe. 

The analysis of our work can be improved in several aspects. In particular, the kinetic 
equations we used are only approximate. They are dealing with a typical particle in the 
plasma with momentum ~ T. Clearly, the exact kinetic equations (to be derived yet) must 
have integro-differential character. We also did not take into account the mass corrections 
O (^-^^ and O (^). We postpone, however, the study of these effects till the moment 
when the singlet fermions are discovered. 
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A. The numerical method 



The straightforward numerical solution of equations (2.10), (2.11) and ( 2.12| ) causes no 
problems for quite a wide range of the parameters of the i/MSM. However, in some corners 
of the parameter space, where the damping is strong and/or oscillations are very rapid, it 
is faced with a number of difficulties. Due to simultaneous presence of many different time 
scales, related to damping of various quantum numbers and due to oscillations between 
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singlet fermions, the corresponding differential equations belongs to the so-called "stiff" 
type. For example, for typical parameters Tw — 100 GeV, M = 1 GeV AMm = 
GeV and e = 0.1, the oscillation rate ~ 2MAMm ^ 2 • 10^^ GeV is some 7 orders of 
magnitude larger than the rate of the singlet fermion creation, ~ sin f-n— 
GeV. So, to solve kinetic equations numerically we will use a method which is very similar 
to the one developed in the Appendix A of Ref. The basic idea is to determine the 

damping and oscillation rates at any given moment of time and then disentangle the kinetic 
evolution accordingly. This is explained in detail below. 

The system of equations ( ^.10| ), ( 2.11 ) and ( 2.12 ) can be rewritten in the form 



-if = W, (A.l) 

I at 

where ^ is a real vector of dimensionality 11, which includes all unknowns, namely 2x4 real 
functions determining density matrices for singlet fermions, 6p^ and 6p^, and 3 chemical 
potentials for active leptonic asymmetries pa- The 11 x 11 complex (but not Hermitian) 
matrix Ti can be called effective Hamiltonian. The (complex, in general) eigenvalues of 7i 
at some moment of time t give the instant relaxation and oscillation rates in the system. 
Let us look for solutions of equation ( |A.lD in the form: 

m = e{t)E{t)Mt), (A.2) 

where il^o{t) is some vector, and e(i) is a matrix build from the eigenvectors of the matrix 
H: 

ne = ep, (A.3) 

where p is a diagonal matrix, constructed from eigenvalues of Ti and the matrix E(t) is 
defined by: 

E{t) = exp - i ^ dt'p{t')^ . (A.4) 

In this way the rapidly oscillating or exponentially damped behaviour of -(/'(i) is factored 
out, and V'o(^) is slowly varying. Now we will derive a convenient form of equation for 

Mt)- 

Let us denote the eigenvectors of by Cj. Then we have Hei = pjCi and e = (ei, 62, ...). 
It is useful to introduce another set of the eigenvectors fj defined by fjTi = fjqj and 

introduce the matrix / as / = I /2 I . The sets of eigenvalues qj coincide with the set pi 

up to permutation, since they are found from one and the same equation, det(T-L — p) = 
or detiTi — q) = 0. If 7^ qj, the vectors and fj are orthogonal, fjCi = 0. Therefore, 
the reshuffling of the set ft can ensure that the matrices / and e satisfy fe = A, where 
yl is a diagonal matrix. We can choose the normalisation of / and e in such a way that 
A = Identity, but we will keep it arbitrary in what follows. 

Using the orthogonality of e and / one can get an equation for ipo{t): 

t - -(-)-'/! -^0, (A.5) 
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and the time evolution for the eigenvectors e{t) and the eigenvalues p{t), f{t): 



% = ^f/^e)., (A.6) 



dt An \ dt 

l^eG, (A.7) 

f^=-AG f, (A.8) 

Gii = 0, Gii = 7 ^-^^(f^e) . (A.9) 

' {pJ-p^)AiiV dt Jij ^ ^ 



The equation (A. 5) is still not suitable for numerics due to explicit presence of expo- 
nentially small terms in the matrix E{t), coming from imaginary parts of the eigenvalues 
and corresponding to strong damping (entering to thermal equilibrium of the corresponding 
processes). To avoid this kind of problem, we remove the imaginary part of the eigenvalues 
from the exponential by the change of variables: 

Mt) = exp ( - ^ dt' Im(p(t'))) (A.IO) 

E{t) = exp - i J dt' Re(p(t'))) , (A.ll) 
V'(t) = e{t)E{t)i;{t) . (A.12) 
The new equation for i/^(t) is: 

^ = Im(p)V; - {AEy'f^Ei;. (A.13) 

These equations were used for compute the behaviour of the density matrices and of the 
chemical potentials with the Wolfram Mathematica command NDSolve. It is convenient 
to use a dimensionless variable z, related to temperature as 



_ hp sin (pMp _ hp sin (pMp 

with hp = 2 X 10^^^ and changing from (initial state) to ~ 1 (sphaleron freezing). 

For some choice of the parameters the computation crashes and further refinements of 
the method were necessary. 

First, if two of eigenvalues are equal to each other at some moment tq, the equation 
( |A.9| ) gets singular. To go around this problem, we continued the coordinate z to the 
complex plane, to go around the singularity: z = t + i ^^^^a^^^"^ ^^-^y where r, a and b are 
some real parameters. The result must be independent of the specific choice of a and b, 
which we verified numerically. 

Second, if the frequency of oscillations between two sterile neutrinos is very high, the 
numerical solution of ( A.13| ) is extremely slow, as a very small time step size is required. 



At the same time, no asymmetry is produced in this regime, exactly due to the fact that 
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the oscillations are rapid 16, 27]. Therefore, we introduced an artificial damping 
of oscillations, effectively averaging them to zero. This was done by multiplying E by 

exp I — A(^J^dt'Re{p{t'))^ ), where A is some real parameter. We verified numerically 
that this procedure indeed works. 



B. Baryon asymmetry as a function of z/MSM parameters 



In this appendix we present some numerical results which clarify the dependence of the 
baryon asymmetry on the parameters of the model, and consider its time evolution. A 
comparison with a perturbative computation is also made. 

For computations, we set the parameters related to the active neutrinos mass matrix 
as specified in eqns. (|33J3]|, |3j) and Tw = 100 GeV. 

In Fig. |5| we present the dependence of asymmetry on the mass M for some choices of 
other parameters (e and AMm)- 




0.5 1.0 2.0 

M [GcV| 



Figure 5: Baryon asymmetry as a function of Af for e = 5 • 10 AMm = 10 ^'^ GeV in blue - 
solid line, for e = 0.5, AMm = 10"^° GeV in red - dashed line and for e = 7 • lO^^, AMm = 10"^ 
GeV in black - dotted line 

In Fig. we show the dependence of the asymmetry on AMm, keeping M = 2 GeV 
and e = 0.5, with latter two numbers chosen ad hoc. The upper line corresponds to the 
perturbative solution given by (|2.16| ), and the lower one to numerical integration. These 

Im> 10"^^ GeV. For smaller AMm 



two curves are almost the same for AMm > 10 GeV. For smaller AMm the conditions 



(ii) and (iii) formulated in Section 2.3 are not valid any more. The asymmetry decreases 
for very small AMm, as the resonance has no time to develop. 

In Fig. |b we fix Af = 2 GeV and AMm = W'^ GeV and vary e. The analytical 
perturbative solution is close to the numerical one for e > 0.01. This is what is expected: 



the condition (i) of Sect. 2^ requires e > 0.015 to insure that the sterile neutrinos are out 
of thermal equilibrium. For smaller e the asymmetry decreases because the singlet fermions 
start to equilibrate. 

In Fig. |6|c we fix e = 0.5, AMm = 10~^ GeV and vary M. The curves start to deviate 
from each other at large masses, where the condition (i) is not valid any longer. 

In Fig. we show the time evolution of different asymmetries for the following choice of 
parameters (for the normal hierarchy of neutrino masses). We take e = 0.015, AMm = 10~® 
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GeV and M = 1 GeV, and consider evolution of the system till temperature T = 100 GeV. 
These parameters are chosen in such a way that equilibration of singlet fermions occurs at 
T > T]y (so that Tjytw > !)• Therefore, one would expect that the asymmetry start to 
decrease after t ~ l/Tj^. Also, for this choice Tl > T\y (see Sect. 2^), indicating that the 
production of asymmetry in A''2,3 occurs around the time corresponding to T^. The time 
variable r used in these plots is defined in appendix 

The asymmetries start to be produced around r ~ 0.05 (see Figs. |^a,b,c and below). 
The asymmetries in singlet fermion is almost the same as in (Fig. ^). The asym- 
metries in /i and r flavours are very close to each other (upper line in Fig. |^b). They start 
to decrease for r > 0.5, which is the onset of damping (see below). The graph (0c) shows 
us that the asymmetry in the active sector is the same with the opposite sign as the one 
in the singlet sector. 

In Fig. 1^ we show the time evolution of the damping (left panel) and ^^^^2 — -^3 oscillation 
factors (right panel) , appearing as exponentials in ( |A.10D , for the relevant eigenvalue of the 
effective Hamiltonian %: Jlm(pE2)dt and J l!ie{pE2)dt- The numerical and analytical 
results of the first expression almost coincide for this particular choice of parameters, but 
not the real part due to our choice of CP-violating phases. 

The oscillation exponent Fig. ^ (right panel) reaches 1 at r ~ 0.05, which is the time 
when the asymmetries are produced on the graphs of Figure 0. The damping exponent 
Fig. ^ (left panel) approaches 1 at r ~ 0.6, corresponding roughly to maximum of the 
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T 

(c) 

Figure 7: Asymmetries in the sterile sector (a) (blue - solid line for N2 and red - dashed line for 
A^s), in the active sector (b) (red - long line for e, blue - short line for fi and brown - dotted line for 
t) and total asymmetries in both sector (c) (blue - solid line for the sterile sector and red - dashed 
line for the active sector). 



asymmetries in Fig. |^. 

In general, the results of numerical integration are in accordance with the qualitative 



picture developed in [27|. 
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Figure 8: Time dependence of J Im{pE2)dt (left) and of J Re{pE2)dt (right). The numerical and 
analytical results are the same (blue - solid line for the numerical results and red - dotted line for 
the analytical solution. 
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C. CP phases and maximal baryon asymmetry 



Our procedure to find the parameter space of the i^MSM leading to acceptable baryon 
asymmetry of the Universe implies the maximisation of the result with respect to three 
unknown phases r], (p and a for normal hierarchy and ry, (p and for inverted. With the 
choice of active neutrino mass parameters ( tj.l|j3.2| ) the numerical analysis shows that it is 



achieved close to r] 



a 



and ^ 



(Note that for 6*13 = and 6^23 



and 



parametrisation of active neutrino mixing matrix as in |28|] the Dirac phase cp plays no role 
for inverted hierarchy.) 

In Fig. 1^ we show the dependence of asymmetry on cj) for a typical choice of M, AMm, 
and e (left panel, fixing r] = a = ^) and on a (right panel, fixing r] = (p = In Fig. 10 
we show the dependence of asymmetry on ^ (fixing r] = a 
hierarchy. 



^) for the case of inverted 





-0.00001 
-0.00002 

-0.00003 



Figure 9: Lepton asymmetry as a function of (f> (left panel) and a (right panel) for normal hierarchy. 




Figure 10: Lepton asymmetry as a function of ^ for inverted hierarchy. 

The boundaries of the parameter space given in Sect. ^ correspond to the extremal 
choice of phases, found in this Appendix and to the observed baryon asymmetry. Below 
we will determine the values of M, AMm and e which extremise asymmetry and find its 
maximum (using, again, the set of phases (|3.3|) ). 

As a guiding line for numerics, we can proceed first with the analytical estimate. In 
Fig. 11 we plot the dependence of 5cp, defined in ( p. 17 ) on e. In general, CP-violation is 
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larger in inverted hierarchy case and is maximal around e ~ 0.5. Replacing the brackets in 
(2.16) by 1, we arrive to 

~ 1.5 X 10~^ (2.5x10"'^) (C.l) 



riB 
s 



for inverted (normal) hierarchy. The asymmetry is expected to be maximal for 

""^-^^^^"^1, 1. (C.2) 

8eTw AMAMmMo ^ ' 




Figure 11: Sep in function of e for the normal hierarchy (blue - solid line) and for the inverted 
hierarchy (red - dashed line). 



These expectations can be verified by numerics. In Fig. 12 we show the dependence 
of the maximal asymmetry on the singlet fermion mass (for each choice of M and Tw we 
find the values of e and AMm that give the maximal asymmetry). 
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Figure 12: Maximal asymmetry as a function of M for Tw = 100 GeV (blue - short dashed line), 
T\Y = 150 GeV (red - dashed line) and Tw = 200 GeV (green - long dashed line). Left panel - 
normal hierarchy, right panel - inverted hierarchy. 

In Fig. |l^ we show for which e and AMm the asymmetry is maximised for different 
T\Y, with the parameter along the curves being the mass M. 



In Fig. 12, the maximal asymmetry for the normal hierarchy is reached for masses 



M bigger than 50 GeV. For the inverted hierarchy, the numerical computation of the 
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Figure 13: The parameters e and AMm corresponding to the maximal asymmetry for Tw = 100 
GeV (blue - short dashed line), Tw — 150 GeV (red - dashed hne) and Tw — 200 GeV (green - long 
dashed line). The mass M grows to the right Left panel - normal hierarchy, right panel - inverted 
hierarchy. 



asymmetry maximum is 10^ times smaller than the estimate ( p.l| ). This difference comes 
from the fact that one can not extremise both expressions in ( |C.2 ) at the same time for 
our choice of temperature. Let us take the case T = 100 GeV for example. The maximum 
(for the inverted hierarchy) is reached for M = 25 GeV, e = 0.5 and AMm = 5 • 10~^^ 
GeV which extremises the first expression of (|C.2| ), but the second is of O(10~^). Looking 
at the analytical expression, we should decrease AMm and get more asymmetry, but we 
are already at the limit of validity of the perturbative expansion because of condition 



of Section 2.3 
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